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Abstract 

The quantum instability of the mean-field theory for identical bosons is shown to be described 
by an appropriate Bogoliubov transformation. A connection between the quantum and classical 
linear stability theories is indicated. It is argued that the instability rate in a system of identical 
bosons must be strongly affected by the nonlinear terms (interactions). In the case of the repulsive 
interactions or strong attractive interactions the instability rate is suppressed. On the other hand, 
a weak attraction significantly enhances the instability rate. The results can have applications in 
the field of Bose-Einstein condensates of dilute quantum gases. 
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Under some conditions, a quadratic Hamiltonian describing a system of identical bosons 
can be diagonalized by the Bogoliubov transformation 1]. The physical sense of the di- 
agonalization lies in a redefinition of the vacuum state and introduction of quasiparticles 
satisfying the boson commutation relations. The diagonalization is used, for instance, in 

n 

description of the excitations in a Bose-Einstein condensate of dilute quantum gas |2J. 

Diagonalization is not always possible which means that the quasiparticle excitations may 
give rise to an instability jj. Below it is shown that such instabilities are described by an 
appropriate Bogoliubov transformation which, however, is not related to the diagonalization 
procedure. Moreover, it is argued that even weak interactions (i.e. the higher-order terms) 
can cause a dramatic effect on the instability rate. 

Driving a quantum many-body system through an instability region by changing the 
svstem parameters is an important tool in the quantum state engineering in BEC physics 
2], for instance, in the process of the atomic association 0,0,0] and molecular dissociation 
in BEC . In a region of instability the mean-field theory breaks down 
instability in system of identical bosons is also referred to as the curve crossing 
related to the Landau- Zener theory . In a condensate with the attractive interactions 
the quantum corrections are crucial for description of the dynamics when the mean field 




state is unstable 

Let us first briefly recall the basics of the diagonalization procedure Consider a 
Hamiltonian quadratic in the boson creation (&j) and annihilation (bj) operators: 

H = Y1 A *Abi + lYl B *M + \ E B *kMi> (!) 

k . / k ; Z k . / 

where A* = A and B = B (here t, ~, and * denote the Hermitian conjugation, the transpo- 
sition, and the complex conjugation, respectively). The diagonalization of the Hamiltonian 
((TJ) introduces new boson operators, say Q and such that the corresponding matrix A 
becomes diagonal and B vanishes. 

The fact that the new operators satisfy the boson commutation relations together with 
the linearity and homogeneity of the transformation gives immediately the most general 
form of the Bogoliubov transformation: 





rfU - V^V = J, UV-VU = 0. 



(2) 
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Here, for simplicity, the index-free operators denote the whole column of such operators, 
e.g. £ = (£1, . . . , £ n ). Note that the conditions on matrices £/ and V in equation (J2J) also 
guarantee the existence of the inverse transformation 

(3) 

To diagonalize the Hamiltonian (0) the Bogoliubov transformation must be such that 
columns of the transformation in equation (JHJ) are the eigenvectors of the following matrix 
eigenvalue problem 



A, ' , (4) 








where U.^ and V.^ denote the k-th columns of the respective matrices. If the diagonalization 
is possible, the Hamiltonian (JTJ) in the quasiparticle representation reads H — ^ AfcrfA- — 

n 

Afc (Yli \ v i,k\ 2 )> w here vi k are the elements of the matrix V Note that the number of 
quasiparticles N = is conserved. 

The diagonalization procedure described above is closely related to the classical theory 
of stability of the stationary points of the Hamiltonian systems Indeed, the second 

matrix on the l.h.s. of equation (jlj), 



(5) 



is the Hessian of the classical Hamiltonian system corresponding to the quantum Hamil- 
tonian (P) when the boson operators b and are replaced by the c-numbers. Hence, the 
eigenvalues A& of the eigenvalue problem (jH) come in quartets {A, A*, —A, —A*}. For the 
diagonalization to exist it is necessary that there are n non-zero real eigenvalues and n 
corresponding eigenvectors with the positive normalization, i.e. WU — V'V = I. 

The simplest sufficient condition for the diagonalization is the positivity of the corre- 
sponding Hessian. In the classical stability theory the quadratic form of the Hessian is the 
first non-zero term in the energy expansion, thus its positivity implies linear stability jl^ . 
The physical sense of the positivity of the Hessian is that the linearized system about the 
stationary solution is equivalent to a system of uncoupled oscillators (in the stationary case). 
In the quantum case the corresponding Hamiltonian, quadratic in boson operators, can be 



diagonalized. However, it is clear that the diagonalization is not always possible (for exam- 
ple, the Hamiltoian (JHJ) below cannot be diagonalized when \S\ < 7). In this case there no 
interpretation in terms of the uncoupled oscillators since the linearized system is hyperbolic. 

A mean-field theory for bosons can be considered stable if the quasiparticle excitations 
remain bounded (for instance, when the Hessian is positive). There is a correspondence 
between the classical and quantum stability. Physically such a correspondence seems to 
be quite clear: due to the linearity of the Heisenberg evolution equations (for b and b') 
the stability in the classical sense implies that the diagonalization of the corresponding 
Hamiltonian, quadratic in boson operators, is possible. In the classical theory it is established 
that under a variation of the parameters the solution to a Hamiltonian system may become 
unstable due to collision of the imaginary eigenvalues with the energy of op pos ite signs (the 
opposite Krein signatures) or due to a collision of the eigenvalues at zero jl2|. Hence, the 
simplest quantum instability is due to a collision of two frequencies at a resonance. 

The simplest example of the quantum instability (or non-diagonalizable boson Hamilto- 
nian) has already appeared in the description of a boson system driven through a resonance 
[3 • The corresponding Hamiltonian reads 0, 0] 

H = [uj + 5(t)]a\ai + [— uj + 5(t)]c4a 2 + 7(0^2 + a i a 2)- (6) 

This Hamiltonian conserves the difference of the number of bosons: [H, a\a\ — a\a2\ = 0. 
The total number of bosons is not conserved which is a reflection of the fact that the 
Hamiltonian © appears as a linearization of the full nonlinear Hamiltonian about the mean- 
field state or in the case of a macroscopically populated source of bosons. For instance, 
such a Hamiltonian appears in the description of the formation of atomic pairs by the 
dissociation of a molecular condensate 0, Q] . In the case of one highly populated molecular 
mode, while the atomic mode populations are not large, the model Hamiltonian is H am = 
[uj + 5{t)]ip\il)i + [—u + 5{t)]ip\il)2 + {g^mfpiiA + fl^m^i^)- ^ a macroscopic number of 
molecules does not change considerably, the annihilation operator for the molecular state can 
be approximated by an order parameter, \l/ m — > (^> m ). Defining 7 = |g(\l/ m )| and i/jj = e lx aj 
with x — ar g{5'(^m)}/2 we obtain the Hamiltonian 
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For the Hamiltonian (JBJ) the corresponding Hessian reads 

/ S + tu 

A B X 

B* A* 



H 



\ 



(7) 



7 

5-u 7 
7 5 + uj 
7 8-u ) 

The characteristic equation is bi-quadratic: A 4 — 2{oj 2 + 8 2 — 7 2 )A 2 + (■y 2 — 5 2 + cj 2 ) 2 = 



V 



it has the following four roots A = ±(u; ± a/5 2 — 7 2 ) (the two ± signs can be arbitrarily 
selected). Therefore the Hamiltonian © can be diagonalized for \5\ > 7 by solving the 
eigenvalue problem Setting a = 5/7 one obtains the diagonalizing transformation 



U 



u x e iei 







u±e 



( —1/2 ^ 

with Mi = (l — [Vo -2 — 1 — o"] 2 ) and m 2 = [Vc 2 — 1 — (without loss of generality, 
the arbitrary phases 0^2 can be set to zero). In particular, the new boson operators read 



■W 2 



V 



M 2 e 




iG 2 



u 2 e Wl 



(8) 



ii = e tdl {uiai - u 2 al), 6 = e ^ 2 ( M i a 2 - ^a{) 



and the Hamiltonian reduces to 



7^ 



5. 



(9) 



(10) 



When cr = 5/7 is time-independent the diagonalization decouples the Heisenberg evolu- 
tion equations for the boson operators £1 and £ 2 : 



d£ 2 



w)6 



(ii) 



Otherwise, there is no point in diagonalizing the Hamiltonian since the Bogoliubov trans- 
formation will be time-dependent and the evolution of the new boson operators £1 and £2 
will be coupled (similar as the evolution equations for the original operators a% and a 2 ). 

In the classical stability theory the Hessian (JJJ) describes a collision of two eigenfrequencies 
at \8\ = 7 (8 — > —7 from below or 5 — > 7 from above) with the appearance of complex 
eigenvalues (eigenfrequencies). This corresponds to the fact that the boson Hamiltonian (JBJ) 
cannot be diagonalized for \8\ < 7. The corresponding quantum instability can be thought 
as the exponential growth of the number of bosons, described by the boson operators a% and 
a 2 , since the total number is not conserved. 



Though the Hamiltonian cannot be diagonalized in the region of instability, the Bogoli- 
ubov transformation (J2J) plays an essential role in the solution of the Heisenberg evolution 
equations. This statement is valid for any quadratic Hamiltonian (JTJ and follows from the 
fact that the action of the corresponding time evolution operator on the boson creation and 
annihilation operators can be represented by an appropriate time-dependent Bogoliubov 
transformation (which is different from the diagonalizing transformation). The latter fact 
is just a consequence of the boson commutation relations for the creation and annihilation 
operators. 

Let us consider the Hamiltonian (jHJ) as an illustrative example. First one notes that uo can 
be eliminated by the unitary transformation with the generator CL]CL\ — CL2CL2 5 

which results in 

the substitution a\ — > e lujt a\ and a 2 — > e~ lujt a 2 . Then, taking into account the symmetry of 
the transposition a\ <-> a 2 , the Bogoliubov transformation solving the Heisenberg evolution 
equations can be cast as 

("} = ("--")(*). uJ'"°), vJ ° -"'A. (12) 
V a 7 \~ V U ) \al) \0 u\) \-u* 2 J 

The condition U'U — V'V = I requires that \ui\ 2 — \u 2 \ 2 = 1, while the second condition 
from equation (J2J) is satisfied identically. Differentiating equation (j!2j) with respect to time 
and expressing the constant operators a and aj, through the time-dependent operators a 
and by using the inverse Bogoliubov transformation one arrives at the following linear 
system for the parameters of the transformation 

= -S(t)) (n*)' (13) 

A solution of system (fH?|) determines the corresponding solution to the Heisenberg equations 
for the boson operators a and a) . In fact, the latter equations can be obtained by the formal 
replacement: U\ — > a\ and u 2 — >• a 2 . Note that the transformation defined by the solution 
of equation ()13|) has nothing to do with the diagonalization of the Hamiltonian. 

System (|T3*|) formally coincides with the normalized classical system, which obtains from 
the Hamiltonian (jH)) by replacing the boson operators with c-numbers (equations ()22)1 and 
(f2*3*|) with g — 0). The Hamiltonian (JHl) can also be cast as a linear combination of generators 
of the group SU(1, 1) (in our case: /C = {a\ai + a 2 a 2 + l)/2, /C + = a\a 2 and /C_ = a\a 2 
with [JC ,1C±] = ±/C± and [/C+,/C_] = — 2/C ), thus the quantum instability case is formally 
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similar to the problem of parametrically forced quantum oscillator |13j . However, the physics 
is quite different. For instance, in the quantum instability case there no adiabatic frequencies 
in the region of instability. The time-dependent Bogoliubov transformation was first used 
in the solution of the forced oscillator problem in Ref. JlJ . 

System (jT5j) is equivalent to the approach of Refs. 3j, |9(. Let us show, for instance, how 
the asymptotic formula of Ref. Q], relating the boson operators at t — > ±00 in the case 
S(t) = xH, appears in our approach. The first order system (j!3|) is equivalent to a second 
order equation (satisfied by both u\ and w 2 ). Setting t = z/a with a 2 = 2ix 2 we get the 
equation for the parabolic cylinder (Weber) functions: jl^ 



d 2 u k {z/a) 
dz 2 



+ 



z 2 11 







(14) 



with v = z7 2 /2x 2 . Using the simplified version of the integrating Bogoliubov transformation 



a 1 




(15) 



the asymptotic conditions U\ — > 1 and U2 — > as ze~ m ^ 4 — > —00 complemented with the 
asymptotics of the first derivative derived from equation ([T3|) and the standard asymptotic 
formulae for the parabolic cylinder functions one can arrive at the asymptotic formula of 
Ref. 0: 

a k (+oo) = e^ 2/2 ^a k (-oo) + e^ 2 /* 2 _ i]i/2 a t / (_ oo)j 



2~ i7 2 /2>t 2 g— i7r/4 



r(i - i 7 2 /2x 2 
r(i + z 7 2 /2x 2 )_ 



1/2 



(16) 



with k' — 2 for = 1 and k' = 1 for k = 2. 

I — I 

It is noted that formula (J16|) resembles the well-known Landau-Zener formula [l6j (where 
the Weber functions also play the key role, see for instance Ref. ^3|). However, the quantum 
instability due to a collision of two eigenfrequencies is completely different from the Landau- 
Zener problem. This is manifested by our system (|13p. Indeed, system (113)1 resembles 
the Landau-Zener system but the evolution matrix on the r.h.s. of equation (jl3|) is not 
Hermitian and the normalization condition for the "wave vector" (2x1,1*2) is formulated in 
the Minkowsky space rather than in the Euclidian one. One consequence is that, as distinct 
from the Landau-Zener case, there is no equivalent of the adiabatic transition for slow 
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time dependence of a = 5/7, since the boson system passes through a region (\8\ < 7) 
where the diagonalized Hamiltonian (fTUj) is invalid. In other words, there is no adiabatic 
eigenfrequencies slowly depending on time in the region of \5\ < 7, instead the system is 



Therefore, the principal physical difference from the Landau-Zener problem lies in the fact 
that, due to the very instability that system ()13|) describes, the number of bosons grows in 
time and their interaction cannot be neglected at all times, in general. Let us consider what 
effect the interactions may have on the instability rate. Remembering that the Hamiltonian 
(JBJ) appears as an approximation to the full nonlinear Hamiltonian (for instance, with the 
binary atomic collision term in the case of molecular dissociation) and eliminating the u 
by cli — > aie tujt and a 2 — > a 2 e~ lujt one arrives at the simplest non-quadratic Hamiltonian 
describing the instability and interactions 

H = 8(t)(a\ai + a\a 2 ) + l(a\a\ + a\a 2 ) + gn{a[ai) 2 + gi 2 a\aia\a 2 + g22{a 2 a 2 ) 2 . (17) 

Below, for simplicity, we will set all the interaction coefficients equal, g = gu — g\ 2 = 
#22- The Heisenberg evolution equations written for the operators di = e~ igKt ai and a\ = 
e~ i9Kt a 2 , with K = a[ai — a\a 2 , read 



where fi\ = a[ai and h 2 = a 2 a 2 . 

System (j!8)) is equivalent to an infinite-dimensional linear system of equations arising in 
the Fock representation of the Schrodinger equation with the Hamiltonian (|17jl . It turns out 
that a satisfactory numerical simulation of such a system requires a large number of the Fock 
modes to be used (about 10 000) due to the spreading of the excitations over a large number 
of the Fock amplitudes. This is in contrast to the linear case (g = 0) when the corresponding 
infinite dimensional system in the Fock space is equivalent to the finite-dimensional system 
flllSjl for the parameters of the the Bogoliubov transformation. 

Therefore, in the nonlinear case one has to rely on an approximate system. Our approx- 
imation is based on the replacement of the number operators hi and h 2 in equation (fTSj) by 
their averages. Then the evolution of the boson operators di and a\ must be given by an 
appropriate Bogoliubov transformation. The parameters of the corresponding Bogoliubov 



hyperbolic with the solution given by a combination of the exponents e 




(18) 
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H 


\u' 2 ) 










transformation satisfy the system 

6(t) + g + 3g(n l ) 7 

-7 -8(t) - g - 3g{n 2 ) 

Equation (|15|1 (with all a replaced by a) provides the expression for the boson operators 
through their initial values. Using it one arrives at the following expressions for the averages: 

(ni) = \ui\ 2 (ni ) + |u 2 | 2 ((n 20 ) + 1), (n 2 ) = |wi| 2 (n 20 > + \u 2 \ 2 ((nio) + 1) (20) 

under the assumption that there is no initial pairing of the bosons: (010^20) — 0. Using 
the normalization condition \ui\ 2 — \u 2 \ 2 = 1 and making the phase transformation Ui = 
e 3ts«n. 20 >+i/2)t w / an( ^ ^ 2 _ e 3is((n 2 o>+i/2)t u / we arr j ve a ^ nonlinear system 

8(t)+g/2 + 3g(N + l)\u' 1 \ 2 7 





H 


I "'2 1 





-7 -5(t)-s/2-3s(iVo + l)K| 2 

(21) 

where JV = (n w ) + (n 20 ). 

At this point it is interesting to compare the approximate system (12 lj) with the classical 
system resulting from the Hamiltonian (fTTjl when the boson operators are replaced by the 
c-numbers: 

h c1 = s(t) [cid + c;c 2 \ + i\c{c; + C X C 2 \ + g[(Cl) 2 Cl + C{C X C*C 2 + (C*) 2 C 2 2 ]. (22) 

The corresponding classical system reads 




C[\ 8(t) + 3gK \C[\ 2 7 

C'l) \ -7 -8(t)-3gK \C' 2 \ 2 




(23) 



where d = e iaKot ^7QC[, C* 2 = e iaKot ^K~ Q C'* 2 , with K = |C X | 2 - |C 2 | 2 being an integral of 
motion. Here \C[\ 2 — \C 2 \ 2 = 1. 

Systems (I23|) and (}2T]) are similar. One difference between them lies in the definition 
of the nonlinearity coefficient. In the classical system it is proportional to the integral of 
motion, while in system (|21j) the nonlinearity coefficient is determined by the initial average 
of the total number of the bosons. (Note that for Kq = the correspondence between 
systems (J2~3*j) and (|2T|) cannot be established.) Another difference between the two systems 
lies in the physical sense of the variables. In the classical system C' k are the (normalized) 



amplitudes of the perturbation about the mean-filed state, while in the case of system (|21|) 
Uk are parameters of the Bogoliubov transformation which describes the evolution of the 
boson operators. 
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FIG. 1: Numerical solution of the reduced system (|21|) in the case of repulsive interactions between 
the bosons. Here *c = O.57. The dashed line is the asymptotic value of (hi) = \u2\ 2 for g = 0, 
according to formula ifTBl) . 



Consider now the effect of nonlinearity on the quantum instability. Let us assume that the 
system is driven through the instability linearly, 5 = x 2 t. Consider, for example, evolution 
of the system which was initially in the vacuum state (i.e. N = (h 10 ) + (7120) = 0) and far 
from resonance: —to 3> l/^ 2 - 

We are interested in the average values of the number operators, which, according to 
formula (J20J), in the case of iVo = read (hi) = (h 2 ) = |m 2 | 2 . 

We have found that in the case of the repulsive interactions between the bosons, g > 0, 
the instability rate is significantly reduced and the asymptotic averages are lower by orders of 
magnitude, see figE Note that even a weak repulsion has a dramatic effect on the resulting 
asymptotic number of the bosons. 

In the case of the attractive interactions, g < 0, the effect of a weak attraction consists 
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in enhancing the production of the bosons, as seen from fig. El On the other hand, a strong 
attractive interaction decreases the instability rate. 




FIG. 2: Numerical solution of the reduced system (|21jl in the case of attractive interactions between 
the bosons. Here k = 7. The dashed line is the corresponding asymptotic result for g = 0. We 
have used to = —100. 

In conclusion, it is shown that the quantum instability of the mean-field theory for iden- 
tical bosons can be described by an appropriate Bogoliubov transformation. The relation to 
the instability in the corresponding classical system is established. We argue that a quantum 
instability in a system of identical bosons is strongly affected by the interactions. In the case 
of repulsion or strong attraction between the bosons the instability rate is suppressed. On 
the other hand, a weak attraction significantly enhances the instability rate. These results 
can have applications in the field of Bose-Einstein condensates of quantum gases. 
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